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ABSTRACT 

Although the spectral lines of hydrogen contain valuable information on the physical properties of 
a variety of astrophysical plasmas, including the upper solar chromosphere, relatively little is known 
about their scattering polarization signals whose modification via the Hanle effect may be exploited for 
magnetic field diagnostics. Here we report on a basic theoretical investigation of the linear polarization 
produced by scattering processes and the Hanle effect in Lya, Ly/3 and Ha taking into account 
multilevel radiative transfer effects in an isothermal stellar atmosphere model, the fine-structure of 
the hydrogen levels, as well as the impact of collisions with electrons and protons. The main aim of 
this first paper is to elucidate the key physical mechanisms that control the emergent fractional linear 
polarization in the three lines, as well as its sensitivity to the perturbers density and to the strength 
and structure of micro-structured and deterministic magnetic fields. To this end, we apply an efficient 
radiative transfer code we have developed for performing numerical simulations of the Hanle effect 
in multilevel systems with overlapping line transitions. For low density plasmas such as that of the 
upper solar chromosphere collisional depolarization is caused mainly by coUisional transitions between 
the fine-structure levels of n = 3, so that it is virtually insignificant for Lya but important for Isyfi 
and Ha. We show the impact of the Hanle effect on the three lines taking into account the radiative 
transfer coupling between the different hydrogen line transitions. For example, we demonstrate that 
the linear polarization profile of the Ha line is sensitive to the presence of magnetic field gradients in 
the line core formation region and that in solar-like chromospheres selective absorption of polarization 
components does not play any significant role on the emergent scattering polarization. 
Subject headings: magnetic fields — polarization — radiative transfer — scattering — Sun: chromo- 
sphere 



1. INTRODUCTION 

This paper is the first of a series aimed at elucidating the physics and the diagnostic potential of the linear polarization 
produced by scattering processes in some hydrogen lines formed in stellar atmospheres (i.e.. Ha, Lya and Ly/3), 
taking into account their sensitivity to magnetic fields via the Hanle effect. This is of p articular interest for the 
exploration of magnetic fields in the outer solar atmosphere fe.g.. iTruiillo Bueno et al1l2005l) . bu t also for facilitating 
the u nderstanding of the spectral line polarization observed in other astrophysical systems (e.g.. [Harrington &: Kuhnl 
l2007t l and the development of suitable polarized radiation diagnostics. 

The true physical origin of the scattering line polarization is the presence of atomic level polarization (population 
imbalances and quantum coh erences between the sublevels of dege nerate atomic levels), caused by the absorption 
of anisotropic radiation (e.g.. 'Landi Deg l'Innocenti fc Landolfil l2003 ) . Typically, in weakly magnetized atmospheres 
anisotropic radiation pumping produces atomic level alignment (i.e., the populations of sublevels with different values of 
\M\ are unequal, with M the magnetic quantum number) and the ensuing selective emission and/or selective absorption 
of polarization components give rise to linear polarization in the emergent spectral line radiation. In principle, atomic 
level orientation may also be induced (i.e., an atomic excitation such that the sublevels with M > are differently 
populated from those with Af < 0), but since its production requires very especial circumstances (e.g., a significant 
net circular polarization in the incident radiation) here we assume that the atomic levels are aligned, but not oriented. 

In the atmospheres of the Sun and of other stars the hydrogen lines considered here are optically thick. Therefore, 
in this investigation we take into account radiative transfer effects in given atmospheric models. The numerical solu- 
tion is obtained by solving jointly the statistical equilibrium equations for the multipolar components of the atomic 
density matrix corresponding to each level of total angular momentum j and the Stokes vector transfer equation for 
each of the allowed t ransitions in the multilevel atornic mo del under consideration. This so-called non-LTE prob- 
lem of the 2nd kind (jLandi Degl'Innocenti &: Landolfil l2004f) can be solved through the application of very e fficien t 
iterative methods and accurate formal solvers of the Stokes- vector transfer equation (jTruiillo Buenol Il999l |2003[) , 
as shown by iManso Sainz fc Truiillo Buend ()2003al l2010t) for the Call IR triplet and by IStepan et all (|2007l) and 
IStepan fc Truiillo Buenol ()2010f ) for the Ha line. 
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Fig. 1. — Grotrian diagram of a three n-level hydrogen model atom (not in scale). The components of Lya, Ly/3, and Ha arc indicated 
along with their experimental Einstein A^i coefficients (in s~^). The solid lines correspond to the transitions involving a polarizable upper 
level. Such transitions play a dominant role in the scattering polarization problem of the above mentioned lines. The arrows indicate the 
collisional dipolar transitions, which play a depolarizing role. 

The problem of scattering polarization in hydrogen lines has attracted the a ttention of various researchers before, but 
concerning mainly the diagnostic problem of solar prominences. For example, iLandLDegl'Innocenti et al. (1987) used a 
simplified radiative transfer model to estimate the expected linear polarization in Ha and H/3 in quiescent prominences 
having a non-negligible optical thickness in Ha. The authors assumed that the /-component of the source function 
corresponding to each hydrogen line is constant within the prominence, which was schematized as an "infinitely sharp 
slab" standing vertically over the solar surface. The solution of the radiative transfer equation for each hydrogen line 
transition, assuming that the ensuing source function is constant, allowed them to compute the spectral line intensity 
at each height within the slab and the radiative rates that enter the statistical equilibrium equations for the multipolar 
components of the atomic density matrix. The numerical solution of the resulting linear set of equations allowed them 
to obtain a zero-order estimate of the /, Q and U emissivities and to compute the ensuing emergent Stokes profiles. 

Concerning the Lyman lines, earlier investigations considered always the optically thin case of the solar corona 
obser ved off-the-limb at large distanc es above it, either neglecting the effects of integration along the line of 
sight dBommier &: Sahal-Brechotlll978l) or taking them into account assuming g iven coronal magnetic field models 
(|Fineschi et al.lll993t iRaouafi et al.ll2009t iDerouich et al.ll201(3l : iKhan et al.ll2011[ ). Anoth e r inte resting o ptically-thin 



invest igation is that of ICasini fc Manso Saind (|2006[ ) who expanded the work of ICasinil ()2005| ) and of iFavati et al.l 
()1987( ) by considering the effects of turbulent microscopic electric fields on the scattering polarization of hydrogen 
lines. Their work was based on the quasi-static approximation for describing the interaction of charged particles with 
the radiating hydrogen atoms, which might not be suitable for quantifying the significance of possible polar ization 
signa tures caused by random electric fields in the low density plasma of solar chromospheric structures (e.g.. IStehlel 
fl990l) . 

When the line of sight (LOS) points to the solar disk the observed spectral line radiation in the Lyman lines does not 
originate in the corona, but in the optically-thick plasma of the solar transition region (or in that of a filament, if such 
a plasma structure happens to be located along the LOS). Therefore, the possibility of having measurable scattering 
polarization in the Lyman lines when doing on-disk observations is of great diagnostic interest, because of its sensitivity 
to th e Hanle effect in the s olar t ransition plasma. In their paper on the scientific case for spectropolarimetry from 
space Iryujillo Bueno et aP ()2005f ) argued that on-disk observations in the strongest lines of the Lyman series should 
show scattering polarization signals, even in forward scattering geometry at the solar disk center when in the presence 
of an inclined magnetic field. That tentative conclusion was based on a simplified model, but in a forthcoming letter 
and in paper-2 of this series we will show through radiative transfer calculations in semi-empirical and hydrodynamical 
models of the solar atmosphere that measurable scattering polarization in the Lyman lines is indeed expected (see 
iTruiillo Bueno|[201Cll for a brief overview). 

The radiative transfer problem of scattering polarization and the Hanle effect in hydrogen lines is, however, rather 
complicated and it has not been systematically studied before. Solving this type of radiative transfer problems requires 
finding the self-consistent solution of the atomic level polarization produced by radiatively induced population imbal- 
ances and quantum coherences, taking into account the Hanle effect in multilevel atomic systems with overlapping line 
transitions. To this end, we have developed and applied the very efficient numerical methods outlined in Appendix lAl 
with which we can investigate a variety of interesting radiative transfer problems in solar and stellar physics. In this 
first paper we focus on the case of an exponentially stratified isothermal model atmosphere, an apparently simple 
but very suitable model to elucidate the physics of formation of hydrogen lines in an optically thick medium as the 
extended solar atmosphere. It is interesting to note that the bulk of the solar chromosphere is essentially isothermal 
because t he dissipatio n of mechanical energy there is taken up by latent heat of ionization, and is rapidly lost by 
radiation () Judge! 120061 ). Therefore, for some spectral lines (e.g., the solar Ha line) the isothermal atmosphere model 
is not as unreasonable as one might think at first sight. Moreover, it might also be a reasonable choice for a first 
estimation of the scattering polarization amplitudes of the Lyman lines in filaments embedded in the solar corona. 
However, the isothermal atmosphere model is certainly unsuitable for predicting the linear polarization signals of Lya 
and Ly /3 produced by scattering processes in the solar transition region plasma itself. As mentioned above, the aim 
of this first paper is not to make predictions on the scattering polarization amplitudes of the hydrogen lines, but to 
gain physical insight to facilitate the understanding of a complicated radiative transfer problem of great diagnostic 



Scattering Polarization of Hydrogen Lines 



3 




-0.002 0.002 

Wavelength shift [A] 




-0.0005 0.0005 0.001 

Wavelength shift [A] 



■0.05 

Wavelength shift [, 



Fig. 2. — Displacements of the line components from the line central wavelength for Lya (left), hy 13 (middle)^ and Ho (right). The 
full bars indicate the transitions with a polarizable upper level (playing a major active role in the line polarization) while the empty bars 
correspond to transitions with an unpolarizable upper level. 
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Fig. 3. — The energy levels of n = 2 (left) and n = 3 (right) as a function of the magnetic field strength. Note that the natural width 
of the 2pi/2 arid 2p^/2 levels, F = A^i/2tt, is 100 MHz and the largest natural widths among the n = 3 sublevels are F = 27 MHz (those 
of 3pi/3 and 3p3/2)- In both cases, the natural width is much smaller than separation of the relevant levels. This remains valid unless the 
perturber's density is very large (A^pert ^ 10^'^ cm~^). 

potential. 

At present, all the above-mentioned physical ingredients can only be taken into account within the framework of the 
quantum theory of spectral line polarization (e.g.. .Landi Degl'Innocenti fc Landolfi 2004) . which treats the scattering 
line polarization phenomenon as the temporal succession of Ist-order absorption and re-emission processes, interpreted 
as statistically independent events. This complete redistribution (CRD) approximation is suitable for subordinate lines 
like Ha (see iStepan fc Truiillo Buenol |2010[ ) and for estimating the line-center scattering polarization amplitudes of 
Ly a and Ly /?, but it cannot be used for modeling the wings of the scattering polarization profiles of the same lines 
(e.g.. ISampoorna et al.|[2bl0l and more references therein). 

This paper is organized as follows. After describing in fj^l the atmospheric and atomic model, we discuss in fjjj the 
transfer equation for the Stokes parameters and the statistical equilibrium equations for the elements of the atomic 
density matrix. Scction|4]discusses in great detail our results for the unmagnctized reference case, with emphasis on the 
observational signatures caused by the radiatively induced polarization of the hydrogen levels. Similar discussions can 
be found in ^for the magnetized case, distinguishing between micro-structured and deterministic magnetic fields. Our 
conclusions are summarized in fJSl while three appendices provide information on the numerical method of solution, 
on the role of collisions and on the equations that correspond to the micro-structured magnetic field case. 

2. THE ATMOSPHERIC AND THE ATOMIC MODEL 

We consider an isothermal, one-dimensional, plane-parallel, semi-infinite, exponentially stratified solar model atmo- 
sphere. The geometrical extension of the model atmosphere is 2500 km and its total optical thickness is 4.5 x 10^" at 
the center of the Ly a line, 7.1 x 10^ at the center of Ly /3 and 1.2 x 10^ at the Ha line center. The kinetic temperature 
is T = IC* K. For simplicity, the calculations are carried out using constant values of the electron and proton densities 
(denoted respectively by ric and Hp). Our main aim is to study in depth the non-LTE problem of the formation, in 
this apparently simple stellar atmospheric model, of the linear polarization profiles of the Lya, Ly/3, and Ha lines of 
hydrogen in the absence and in the presence of "weak" magnetic fields (e.g., with magnetic strengths i?<100 G). 

Under such circumstances the linear polarization of Ly a, Ly /?, and Ha is fully dominated by atomic level polarization 
and the Hanle effect, with no significant impact of the transverse Zeeman effect, because the wavelength shifts (produced 
by the Zeeman effect) between the tt {AM = Af„ — AIi = 0) and (7^ ,, {AM = ±1) components {M being the magnetic 
quantum number) are only a very small fraction of the width of the profiles. In this first paper we include only the 
thermal broadening of the spectral lines and we do not consider any background continuum opacity. All the ingredients 
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TABLE 1 

Critical field of the 
hydrogen levels. 



Level Lande factor Bu [G] 

2p3/2 1.33 53.5 

3p3/2 1.33 16.2 

3^3/2 0.80 9.2 

3^5/2 1.20 6.1 



needed for radiative transfer modeling of spectro-polarimetric observations (e.g., the Stark broadening effect) will be 
included in our forthcoming papers. 

Figure [T] shows the hydrogen atomic model we have used in this study, indicating with the usual notation the n£j 
values of each level (with n being the principal quantum number, i the orbital quantum number and j the total angular 
momentum). The model atom includes the 9 fine-structure levels of the first three n- levels of neutral hydrogen, between 
which 11 allowed radiative transitions take place: 2 for Lya, 2 for Ly/3 and 7 for Ha. The energies of these levels 
and the Einstein A^i coefficients of the allowed transitions have been taken from the NIST online database'^. We 
have checked that similar results for the fractional linear polarization of these lines are obtained when adding the fine- 
structure levels of the n = 4 level. Fig. [2] shows the strengths and wavelength positions of the Lya (left panel), Ly/3 
(middle panel) and Ha (right panel) components. Although in laborat ory experiments such individual fine-structure 
components can be resolved using Doppler-free saturation spectroscopy (ISchawlowlll982D . they are however blended in 
stellar atmospheres. The hyp erfine structure of t h e hyd rogen levels is expected t o play a minor r o le on t he scattering 
polarization (e.g., see iBommier fc Sahal-BrechotI (|1982j ) for the case of Lya and iBommier et ahl (|1986al) for the case 
of Ha). 

Table [1] gives the Lande factors of the levels of Fig. [T]that can carry atomic alignment (i.e., those with j > 1/2), 
calculated using the LS coupling approximation. The same table shows also the Hanle critical field {Bh, in gauss) 
of each of such j-levels (i.e., the magnetic strength for which the level's Zeeman splitting equals its natural width) 
without accounting for the quenching effect of collisions. Note that 

B,^hl!pl^. ,1) 

where tufe and gj are, respectively, the lifetime (in seconds) and the Lande factor of the j level under consideration. 
In this paper we neglect quantum coherences between the magnetic sublevels pertaining to different j-levels. For the 
zero-field case this should be a reasonable approximation because the energy separation between such levels is larger 
than their natural width, even when we take into account that the lifetime of the levels is reduced by collisions with 
protons and electrons (hereafter the perturbers) at plasma densities not very much larger than those of the quiet solar 
chromosphere. In the presence of magnetic fields the ensuing level crossings and repulsions may have a significant 
impact on the above-mentioned j-j interferences. Fig. [3] shows the energy splitting of each j-level due to the presence 
of a magnetic field, both for the n = 2 levels (left panel) and for the n = 3 levels (right panel). It is very important to 
point out that the selection rule A£ = ±1 for radiative transitions inhibits couplings between levels n£j and n(£± l)j'. 
Therefore, we only have to worry about the possibility of level crossings and repulsions between levels n£j and n£j' 
and between levels n£j and n{£ ± 2)j'. As seen in the figure, it seems reasonable to neglect them for the magnetic 
strengths considered in this paper (i.e., B < 100 G). 

3. THE EQUATIONS 

The aim of this section is to provide a brief summary on the equations we have solved for the radiative transfer 
investigation of this paper. Information on the numerical methods of solution can be found in Appendix 13 while 
Appendix [C] describes the efficient strategy we apply for dealing with micro-structured magnetic fields when considering 
multilevel systems. 

3.L The multipolar components of the atomic density matrix 

We quantify the ex citation state of each n£j level by means of the multipolar components of the atomic density 
matrix ()OmoMll977[ ). denoted here with the symbol Pqij) (with K and Q integer numbers such that 0<K<2j and 

—K<Q<K). The pg elements with Q = are real numbers given by linear combinations of the populations, pj (M, M), 
of the various Zeeman sublevels n£jM corresponding to the level of total angular momentum j. The total population 
of the atomic level is quantified by y/2j + Ip^, while the population imbalances among the Zeeman sublevels are 
quantified by p^ (e.g., pg(j = |) = kiPjih |) + - Pjihl) ~ Pji~h~l)))- However, the p^ elements 

with Qy^O are complex numbers given by linear combinations of the coherences, pj{M, M'), between Zeeman sublevels 
whose magnetic quantum numbers differ by Q. In fact, since the density operator is Hermitian, we have that for 
each spherical statistical tensor component Pq with Q > 0, there exists another component with Q < given by 



^ |http: //www.nlst . gov/physlab/data/asd. cfm 
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p^Q = (— 1)'^[Pq]*, with denoting the complex conjugate. In quiet regions of the solar atmosphere it is a good 
approximation to assume that the net circular polarization of the incident radiation is zero; therefore, the odd-K pg 
elements (i.e., the orientation components) vanish. Therefore, the 9-level hydrogen model of Fig. [T] requires 38 pg 
elements to fully quantify the atomic excitation state. 

It is important to note that in this paper the Pgij) elements are defined in a reference system whose Z axis (i.e., 
the quantization axis of total angular momentum) is along the stellar radius vector (i.e., the vertical direction in the 
case of a plane-parallel atmosphere). 

3.2. The transfer equation for the Stokes parameters 

In general, the transfer of polariz ed radiation is described by the following vectorial transfer equation (e.g., 
iLandi Degrinnocenti fc Landolfill200l 

^I = e KI, (2) 

as 

where s is the geometrical distance along the ray under consideration, / = (/, Q, U, V)"^ is the wavelength-dependent 
Stokes parameters, e = (e/, eg, ej/, ey)"'" is the emission vector resulting from spontaneous emission events, and K is 
the 4x4 propagation matrix whose coefficients are rjj (the diagonal element), riQ,rii/, rjv (which account for dichroism) 
and Pq, pu, Pv (which describe the anomalous dispersion) . 

The general expressions of the components of the emission vector e and of the propa gation matrix K are very 
involved and they will not be written here^ (see Eqs. 7.15 in ILandi Degl'Innocenti fc La ndolfi 2004). In the Hanle 
effect regime considered in this paper they are given in terms of the local values of the Pg (j) density-matrix elements 
corresponding, respectively, to the upper (j = j„) and lower (j = ji) levels of the line transition under consideration. 
We point out that only the alignment components with K = 2 contribute to eg, eu, r]Q, riu, pq and pu- The ey, 
r]v and pv expressions depend on the orientation components with K — 1, but in our case the odd-K pg elements 
vanish because we assume that there is no net circular polarization in the incident radiation (see i )3.ip . In solar-like 
atmospheres the expressions for e/ and rji are dominated by the pg components (the overall population of the upper 
and lower levels, respectively), but they also depend on the K = 2 components. The frequency dependence of all such 
radiative transfer coefficients is established by the Voigt {(p^) and Faraday- Voigt profiles, whose dependence on 
the magnetic quantum number M we have neglected here arguing that for i?<100 G the Zeeman splittings of the 
hydrogen levels of Fig. [1] are much smaller than the thermal widths of the spectral lines under investigation. This is 
a very good approximation for calculating the linear polarization profiles of Lya, Ly/3 and Ha, because in the quiet 
solar atmosphere the contribution of the transverse Zeeman effect is insignificant compared with that caused by atomic 
level alignment and its modification by the Hanle effect. As mentioned above pg — 0, so that there is no contribution 
of atomic level orientation to ey, rjy and py. Therefore, the only contribution of atomic alignment to the emergent 
circular polarization comes from the (second-order) term —puQ + PqU of the transfer equation for Stokes V, but such 
a contribution is negligible compared with that produced by the longitudinal Zeeman effect which is not considered in 
this paper. Likewise, the second-order terms —pvU -\- puV and —pqV + pvQ of the transfer equations for Q and U , 
respectively, do not play any significant role in the weakly polarized atmosphere. 

In our calculations of the linear polarization of hydrogen l ines in the model atm osphere we always solve the full 
transfer equation ([2]) via the DELOPAR method proposed bv lTruiillo Buenol (|2003f ). but it is useful to note that the 
following Eddington-Barbier approxima te formula can be use d to estimate the emergent fractional linear polarization 
at the center of a strong line transition (jTruuUo Buenol[2003| ) 

7 -51-^(^0], (3) 

where CTq Po/Po must be evaluated at height in the model atmosphere where the line-center optical distance along 
the line of sight (LOS) is unity a nd w)^, is a numerical coefficient th at depends only on the level's angular momentum 
values j and j' (see Table 10.1 in ILandi Degl'Innocenti fc Lando lfi"2004'). 

The first term in the square brackets of Eq. ^ is due to selective emission of polarization components, caused 
by the population imbalances of the upper level, while the second term is due to selective absorption of polarization 
components (or "zero-field dichroism "), caused by the population imbalances of the lower level. It can be shown (e.g., 
iTruiillo Bueno k Manso Sainj|1999l ) that for resonance line transitions in a weakly anisotropic medium like that of the 
solar atmosphere the crQ(j„) values of the above expression can be estimated from the following approximate expression 

4iju)^v/S,^, (4) 



where 



Jq = I dx (f) ^ (pxixn , (5) 



For the case of non-overlapping transitions see IManso Sainz &: Truiillo Buenol 1120101) , who studied the radiative transfer problem of 
scattering polarization and the Hanle effect in the IR triplet of ionized calcium considering multilevel atomic models. 
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is the familiar mean radiation field intensity, J, of the unpolarized non- LTE problem fe.g.. lMihalaill978[) and 

/f 1 
y — [(3/i' - + 3(m2 - l)Q,o] , (6) 

with (j)x being the absorption profile and x the frequency distance from line center measured in units of the Doppler 
width. 

3.3. The statistical equilibrium equations 
The multipolar components of the atomic density matrix are govened by the following rate equations 

d 



Hanlc 



Radiation 



(7) 

Collisions 



where the first two terms on the rhs are the contributions to the rate of change of Pgij) due to the Hanle effect 
and to radiative transitions (see Eq. 7.78 in lLandi Degl'Innocenti Sz Landolfill2004l ). while the last term is that due to 
collisions (see Eq. 7.101 in lLandi Degl'Innocenti fc Landolfill2004f) . In order to calculate such pg unknowns at each 
grid point of the chosen stellar atmosphere model we have assumed statistical equilibrium, i.e., 

d 



^pSO)^0. (8) 

We point out that since the resulting system of equations is not linearly independent one of the equations (e.g., the 
one for the ground-level population) must be substituted by the trace equation of the density matrix, which establishes 
the conservation of the number of particles 

J2V2J7TTp''o{J^) - 1 . (9) 

ji 

Note that the term of Eq. ([7]) due to radiative transitions includes transfer rates due to absorption (Ta), spontaneous 
emission (Te) and stimulated emission (Ts) from other levels, and the relaxation rates due to absorption (Ra), 
spontaneous emission (Re) and stimulated emission ( Rs) towards other levels. The explicit e xpressions for all these 
transfer and relaxation rates can be found in §7. 2. a of iLandi Degl'Innocenti fc Landolfil ()2004l ). Of particular interest 
is the transfer rate Ta due to absorption from lower levels. The expression for Ta is 

TA{ai3i-KiQi ^ aj-KQ) = {2ji + l)B{ai]i ^ aj) ^ (_i)A'.+Q. ^3(2i^ + i)(2X, + 1){2K, + 1) 




K Ki K,\ ,u 



-QQi-Qr- 



Jo:i''cj,»ai)^ (10) 



where 



J|;(i^ai,a,jj = / diy(l)^j{iy~iyo,j,adi)JQ;H, (11) 



with = 0,1,2 and = -K,,...,K,. Note that JqI{v) are the monochromatic radiation field tensors, defined 
as angle averages of the Stokes parameters multiplied by the geometrical irreducible tensors Tq (see Eq. 5.157 of 
ILandi Degl'Innocenti fc Landolfill2004[ ). while Jq' (i^aj.aiji) are the laboratory frame expressions of the radiation field 
tensors, with (f'iji'^ ~ i^aj,aiji) being the absorption profile of the allowed fine-structure transition under consideration 
(see their explicit expressions in lTruiillo Bu eno 2001). As mentioned above, Jg is the familar J-quantity of the standard 
non- LTE problem. The Jq tensor given by Eq. ([6]), which in solar- like atmospheres is dominated by the contribution 
of the Stokes / parameter, quantifies whether the local illumination of the atomic system is preferentially vertical 
(Jq > 0) or horizontal (Jg < 0). The real and imaginary parts of the Jq tensors (with Q — ±1,±2) quantify the 
breaking of the axial symmetry of the spectral line radiation through the complex azimuthal exponentials that appear 
inside the angular integrals. Obviously, a deterministic magnetic field inclined with respect to the symmetry axis of 
the pumping radiation is needed in order for J'^^ and J^2 to be non-zero in a plane-parallel or spherically symmetric 
model atmosphere. 

4. THE UNMAGNETIZED REFERENCE CASE 

Our three n- level hydrogen atomic model has two strong resonant lines (Ly a and Ly /3) and one subordinate line 
(Ha). It is well known that coupling of th ese lines gives r ise to effects t hat cannot be explained using a two- level 
approximation for the individual lines (e.g. I JefferiesI 119681 : iMihalasI I1978D . In particular, coupling of Ha and Ly/3 
through the common upper level n = 3 affects the formation of these lines via degradation of the Ly /3 photons into Ha 
and Lya photons. Furthermore, each of the excited n-levels of hydrogen consists of several close-lying fine-structure 
levels. These levels are connected by permitted optical transitions with only few of the fine-structure levels of different 
n' . Since the fine-structure levels of each n level are coupled together by the An = collisions, the formation of the 
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Fig. 4. — The line source functions normalized to the corresponding Planck function. The perturber's density is Npert = 4 X 10^" 
Thick lines: all relevant collisional processes are taken into account. Thin lines: collisions between the fine-structure levels in n = 
ignored. 

lines depends on the rate at which thes e collisions are able to shuffle the optical electrons among the fine-structure 
levels of any given n (see Appendix IB.2|) . This process leads to photon conversion among the transitions affecting the 
line source functions and, in the limiting case of strong collisions, the n£j and n£'j' levels behave as a single n-level. 
Moreover, these inelastic collisions connecting nearby j levels tend to depolarize the atomic levels and, consequently, 
modify the scattering pola rization of the emergent radiation (cf., i Bommier et al.iri986bt iBommier fc Sahal-BrechotI 
[19911 ISahai -Brechot et al.l [1996) . As pointed out in Appendix [B] elastic collisions between the radiating hydrogen 
atoms themselves do not play any significant depolarizing role for hydrogen in the low density plasma of the solar 
chromosphere. 

In this section, we discuss the role of several phenomena affecting the formation of the hydrogen lines in our non- 
magnetized model atmosphere. Special emphasis is given to the role of collisions with the ambient perturbers. Among 
other effects we discuss, is the fact that the atomic polarization of the n = 3 levels, due to the anisotropy of the 
incident Ha radiation, has a noticeable impact on the Ly /? radiation, so that its anisotropy turns out to be significant 
in the layers where Ly /3 is optically thick. We point out a possible influence of this effect on the linear polarization of 
the Ly/3 wings. 

4.1. Line source functions 

A detailed discussion of the formation of the unpqlarized hydrogen Ly a — Ly /3 — Ha system in an isothermal atmo- 
sphere can be found in the literature (e.g., lJefFerieslll968t lMihalaslll978f ). Here we briefly review the behavior of the 
line source functions. In addition, we consider the fine structure of the levels and we discuss the role of the An = 
collisional transitions on the line source functions. 

Since the vast majority of the Hi atoms is in the ground state, the optical thickness of the atmosphere is much larger 
in the resonance lines than in the subordinate line. As a result. Ha is formed deeper in the model atmosphere and Lya 
and Ly /3 are optically thick in the layers where Ha becomes optically thin. Since the levels 3pi/2 and 8^3/2 are common 
upper levels of Ha and Ly/3, there is an efficient mechanism of conversion of Ly/3 photons into Ha photons. This 
systematic degradation of the Ly f3 photons is responsible for equality of the Ha and Ly (3 source functions (normalized 
to the Planck function) in the atmospheric layers where the branching ratio ^32/(^31 + ^32) is higher than the escape 
probability of the Ly /? photons (where Ann' is the Einstein coefficien t of spontaneo us emission between levels n and 
n'.) Above such layers, the line source functions become uncoupled ()Mihalajll978l ). In conclusion, Ly/3 is not fully 
thermalized in relatively deep layers where its optical depth would suggest it. On the other hand, Lya is thermalized 
in such deep layers because the level n = 2 is not drained by any subordinate line. 

If the polarization of the atomic levels is taken into account, then e/ and 777 become angle- dependent (e.g., Eqs. (7.16) 
of iLandi Degl'Innocenti fc Landolfill2004[ ). However, if the atomic polarization is small (|Pq'*°| <C |po|) this effect is 
of relatively small importance. Neglecting stimulated emission and the wavelength separation of the line components, 
we can obtain an approximate expression for the line intensity source function, 5*/, in terms of the K = multipoles 
of the atomic density matrix. 



Si 



(12) 



where u = n„£„j„ and I = niiiji denote the upper and the lower levels of the line transition, respectively, while A^i 
and Biu are the Einstein coefficients for spontaneous emission and for absorption respectively. 
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Fig. 5. — Profile-averaged fractional anisotropy of the radiation. The solid lines correspond to the 7 individual components of the Hq 
line (cf. Fig. [2]l. The components of hy a and Ly/3 cannot be distinguished at the resolution of the plot. The vertical arrows indicate the 
heights where the line-center optical path is unity. The thin arrows correspond to disk-center observations (line of sight inclination /j = 1) 
while the thick arrows correspond to close to the limb observation {fj. = 0.1). The perturber's density is Alport = 4 X 10^" cm~^. 

The thick lines in Fig. |4] show the self-consistently calculated line source functions of the three lines versus height 
in the model atmosphere. These source functions have been calculated taking into account all the An^O and An = 
collisions at a uniform electron and proton den sity of 4 x 10^° cm~^ . For such a density value, which is typical of 
the upper chromosphere of the quiet Sun (e.g. iFontenla et al.|[l993[ ). collisional depolarization of the 2^3/2 level is 
practically negligible fe.g., iSahal-Brechot et al.iri996[ ). In order to explore possible effects of the An — colhsions on 
the line source functions, it is useful to compare the results with a model in which these collisions are neglected. The 
line source functions of the model with An = collisions neglected in the n = 3 levels are shown by the thin lines of 
Fig-Sl The Ha source function is slightly smaller in the upper layers compared to the depolarized case, while the source 
function of Ly /? is a bit larger. In the model with An — collisions neglected in the n — 3 levels the populations of 
the 3pi/2 and 8^3/2 levels, which are the upper levels of Ly/3, are slightly larger. Likewise, the populations of the 8^3/2 
and 3(i5/2 levels, which play a key role in the Ha line formation, slightly decrease. The explanation is the following. 
Collisions between the fine-structure levels tend to balance the populations of the n-sublevels. As a result, the 3p term 
populated mainly by the Ly /3 radi ation is collisio nally drained in favor of the 3d one (this is the well known mechanism 
of photon conversion discussed bv'Mihalasl ll978D . 

To some extent, the quantitative results of this section depend on the collisional rates of the individual inelastic 
n£j — n'l'j' (with n 7^ n') transitions. For more details about the collisional rates see §4.51 and Appendix IbI 

4.2. Anisotropy of the spectral line radiation 

The atomic level polarization and, consequently, the polarization of the emergent spectral line rad iation depends 
on th e symmetry properties of the radiation field in the atmospheric region of line formation (e.g., iTruiillo Buenol 
|2001| ). In a cylindrically symmetric atmosphere the key quantity is A = Jq/ Jq (see Eqs. [5]and|6]), which quantifies 
the degree of anisotropy of the spectral line radiation under consideration. Typically, A is zero in the deep layers of 
the atmosphere while it is significant around and above the atmospheric height where the line center optical depth is 
unity. The sign and size of A depends on the gra dient of Sj, with ^ > if the illumination is predominantly vertical 
and ^ < if it is predominantly horizontal (e.g., ITruiillo Bueno|[200lt iLandi Degl'Innocenti fc Landolfill200l) . 

Since the hydrogen lines are composed of several overlapping fine-structure transitions, we may expect the anisotropy 
of the transitions pertaining to the same line to be similar. Fig. [5] shows the degree of anisotropy of each individual 
transition as a function of height in the model atmosphere. For each transition the ensuing A value increases sig- 
nificantly above the height of unity optical depth, due to the lack of incoming photons from the outer layers. The 
anisotropy of the resonant doublet that produces the Ly a line remains virtually zero until relatively large heights in 
the model atmosphere (« 1800 km). The anisotropies in both Ly a components are virtually identical because of their 
small wavelength separation and because they share the same lower level (see the left panel of Fig. [2]). The same 
applies to the Ly/3 doublet (middle panel of Fig. [5]). The difference with respect to Lya is that there is a small but 
significant anisotropy in Ly f3 in the intermediate atmospheric layers where the Ly /? optical depth is much larger than 
unity. As seen in Fig. [Sj the Ly /3 anisotropy in the intermediate layers of the atmosphere increases at the same heights 
where the anisotropy of Ha suddenly increases. This behavior will be explained in the following subsection. Fig. [5] 
shows that the differences among the anisotropies of the different components of the Ha line are small but noticeable. 
This is because the separation of the Ha line components is of the order of 0.1 A (see right panel of Fig. [2]) which is 
comparable to the Doppler width of the line. The individual Ha components are therefore active in slightly different 
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Fig. 6. — Fractional alignment of the levels and influence of collisions and lower level polarization. Top: Lower level polarization is taken 
into account. Bottom: Lower level polarization is suppressed. Left: Calculation with a perturbers density Npert = 4 X lO^'^cm"^. Right: 
Depolarizing collisions in n = 3 are neglected. 

parts of the line. 

4.3. Atomic level polarization 

As mentioned in f}3TT]our choice for the quantization axis, Z, of total angular momentum is along the normal direction 
to the plane-parallel atmosphere (hereafter, the atmospheric frame). Given that in this section we are considering the 
unmagnetized case, the ensuing cylindrical symmetry of the atmosphere implies that the multipolar components 
PQ^Q{n£j) = 0, so that the state of any given atomic level can be described by its overall population, pQ{n£j), and 

the atomic alignment components, {n£j). Moreover, K > 2 multipoles of the density matrix (i.e., Pq {n£j) with 
K = 4,6, . . .) of the levels with j > 5/2 can be safely neglected in our weakly anisotropic stellar atmosphere model 
because they are very small compared to the second-rank components. The top left panel of Fig. |B| shows the variation 



3d 



3/2, 



and 3^5/2. In the remaining part of this secti 



with height of the fractional alignments of the levels 2p3/2, 3p3/2, 
on we consider each of such levels separately. 

2p3/2- This is the only polarizable level of the Lya line. Since the optical thickness in Lya is the largest of 
the hydrogen lines, its anisotropy is negligible up to the uppermost layers of the atmosphere (see Fig. The 
fractional alignment of the 2p3/2 level follows the Lya anisotropy as if it were an upper level of a two- level atom (e.g. 
iTruiillo Bueno fc Manso Sainj|T999l ). Because for the chosen density of perturbers (4 x 10^° cm^'^) the depolarizing 
collisions for n = 2 are inefRcient, the surface value of the 2p3/2 fractional alignment approximately corresponds to the 

theoretical prediction of Eq. HJ i.e., Pg/Po ~ 0-05 because w'^^^ — 0-707 and Jq/Jq ~ 0.07. 

3p3/2- This is a common upper level of Ha and Ly/3. The atomic polarization of this level is affected by optical 
pumping in both lines. In the layers where the model atmosphere becomes partially transparent to Ha photons 
(i.e., 1200 km ^h< 1700 km; see Fig. [5]), the atomic polarization of this level is increased due to the absorption of 
anisotropic radiation in the Ha line. The fractional alignment of the 8^3/2 level is smaller than that of the 8^5/2 or 
3d^/2 levels because the 3p3/2 level is strongly populated by the photons of the optically thick Ly/3 line, which are 
mostly trapped or degraded into Ha photons (see the discussion in At heights h > 1700 km the anisotropy of the 

Ly/3 radiation increases (Fig. [S]) and the 8^3/2 level is aligned directly by the ensuing anisotropic Ly/3 pumping. As 
shown in the top left panel of Fig.[6l the surface value of Po(3P3/2)/Po (8^3/2) lies below the theoretical value (« 0.06) 
predicted by the 2-level atom approximation of Eq. |4| This is partly because the level is radiatively coupled to 2si/2, 
in addition to lsi/2, but mainly because the An = depolarizing collisions play a more significant role in the n — 3 
than in n = 2 fine-structure levels. 

3d5/2- This is an upper level of Ha which is not radiatively coupled to the ground level but it is radiatively coupled 
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to the 2p3/2 level by the strongest transition ol the Ha line (see Fig. [T]). The 8^5/2 level becomes significantly polar- 
ized above the atmospheric height where the anisotropy ol the Ha radiation increases (i.e., above /iwllOO km). The 
fractional alignment ol the 3c?5 /2 level is approximately constant in the intermediate layers ol the model atmosphere 
where Ha becomes optically thin but the Lyman lines are still opaque (i.e., between 1200 and 1700km). The polar- 
ization ol the 3(^5/2 level is decreased by its strong collisional coupling with the 8^3/2 level. However, once the Ly/3 
line becomes optically thin and much more anisotropic, its upper level 8^3/2 is significantly more polarized and the 
collisional transfer ol alignment from 8^3/2 to 8^5/2 eventually leads to an increase ol the 8(^5/2 polarization (see the 
lower left panel olFig. [5]). Since the 8^5/2 level is strongly radiatively coupled to the 2^3/2 level, its polarization is, in 
the uppermost layers of the atmosphere, affected by the atomic alignment ol the 2^3/2 level. Collisional coupling and 
the role ol lower-level polarization in the Ha line will be discussed in greater detail below. 

3d3/2- This level is radiatively coupled to the 2pi/2 level and, via a weaker transition, to 2p3/2- The former transition 
makes the major contribution to the Ha line in the blue part ol the profile while the later contributes slightly in the red 
part (see the right panel ol Fig. [2]). Note that even though the strongest radiative transition involving the 8^3/2 level is 
(2) 

2pi/2~3(i3/2, with 1^3/2 1/2 ~ 0.707, the fractional polarization ol the 8^3/2 level throughout the atmosphere is smaller 

than that ol 8^5/2, whose radiative excitation is solely due to the 2^3/2-8^5/2 transition whose w^f2 3/2 ~ 0.529. This 

is because ol the 2^3/2-8^3/2 transition with the negative lactor ^3^2 3/2 — —0.566 leading to an overall decrease ol 
the upper-level alignment. For a similar reason the alignment ol the 8(^3/2 level decreases in the uppermost layers ol 
the atmosphere where Po(3P3/2)/Po('^P3/2) I'ises. The weak collisional transition 87)3/2-8^3/2 has a similar depolarizing 
effect on 8c?3/2.^ 

The non-negligible amount ol anisotropy in Ly/3 between 1200 and 1700km, where the Ly/3 optical depth is very 
large, can be explained by the above-mentioned coupling ol the Ly/3 and Ha lines. The 87)3/2 level is pumped by 
both Ly/3 and Ha radiation and it is collisionally coupled to the 8CZ5/2 level. The Ha radiation is no longer isotropic 
at heights above 1200 km and the upper levels ol Ha, including 87)3/2, are aligned. The larger the Ha anisotropy the 
greater the 87)3/2 fractional alignment. Since the Ly/3 photons are emitted from the polarized 87)3/2 level, the emitted 
Ly (3 radiation there is sufficiently anisotropic and polarized so as to give rise to a significant contribution to its Jg 
radiation tensor. 

The radiation anisotropy in the Ly /3 line in such relatively deep layers ol the atmosphere might have an impact on the 
polarization ol the Ly /3 line wings. However, in this paper our calculations are restricted to the CRD approximation. 
Within this approximation, the Ly /3 anisotropy in such deep layers implies a polarization signal in the line wings 
(cl. i i4.6|) . Although we cannot make any reliable predictions concerning the scattering polarization in the Ly/3 wings 
where the CRD approximation is unsuitable, the mechanism described in the previous paragraph may be important 
lor luture PRD theories ol multi-level scattering polarization. 

4.4. Influence of lower-level polarization on the Ha line 

II the lower level ol a line transition is polarized then dichroic effects (i.e., differential absorption ol polarization 
compone nts) may affect the polarization ol the emergent radiation, even in the absence ol any sig nificant Zccman 
splitting ()Truiillo Bueno fc Landi Degrinnocentilll997l:muiillo Buenolll999l:lTruiillo Bueno et al.H2002l) . An interesting 
example ol this "zero-field" dichroism can be found in the scattering polarization ol the infrared triplet ol Can 
(iManso Sainz fc Truiillo Buenoll2003bf) and in many other spectral lines ol the second solar spectrum (jTruiillo Buenol 
I2009D . Since the ground level ol Hi cannot be aligned when the hyperfine structure ol hydrogen is neglected, all the 
lines ol the Lyman series are non-dichroic. On the other hand, the linear polarization ol the subordinate lines can, 
in principle, be modified by selective absorption ol polarization components because such lines have always a lower 
fine-structure level with j; > 3/2. In the atomic model ol Fig. [l]the only lower level that can in principle be aligned 
is the 27)3/2 level, which is one ol the three lower levels ol the Ha line. Therefore, Ha is the only spectral line ol the 
atomic model ol Fig. [1] whose emergent linear polarization could be influenced by selective absorption ol polarization 
components. 

The dichroic transitions ol Ha are 27)3/2-8s]^/2, 27)3/2-8(^3/2, and 27)3/2-8^5/2. In our isothermal model atmosphere, 
the polarization ol the 27)3/2 level is only significant in the uppermost layers, where the Lya radiation is anisotropic 
(i.e., above « 1800km, as shown in Fig. [5]). Since Ha is optically thin in such outer layers ol our model atmosphere 
its emergent Q / 1 profile is not influenced by "zero- field" dichroism. 

It is important to note that, in addition to the possibility ol having "zero-field" dichroism, the presence ol lower-level 
polarization may pro duce a significant feedback on the atomic polarization ol the other levels (e.g., Eqs. 81 and 82 ol 
iTrujillo Buenol [200 It ). In order to investigate the impact ol the presence ol atomic alignment in the 27)3/2 level on the 
atomic polarization ol the upper levels ol the Ha line, we have solved the very same non-LTE problem ol the 2nd kind 
lor the atomic model ol Fig. 1 but forcing the lower level 27)3/2 to be completely unpolarized throughout the model 
atmosphere. The resulting fractional alignment ol the hydrogen levels is shown in the bottom left panel ol Fig. [6l The 
infiuence ol lower-level polarization can be valorated by comparing the bottom left panel ol Fig. [6] with the top left 
panel ol the same figure in which the lower-level polarization is taken into account. 

^ Note, however, that concerning collisional transitions we deal with the transfer of atomic alignment by isotropic collisions while 
con cerning radiative transition s we deal with absorption of anisotropic radiation from levels that do not necessarily have to be polarized 
fsee lSahal-Brechot et aLlll996l . for a discussion of the alignment transfer mechanisms in the hydrogen fine-structure levels). 
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Fig. 7. — The normalized line intensity (left) and fractional linear polarization (right) profiles for a close to the limb LOS {fj, = 0.1) 
calculated for several densities of the charged particles. The individual panels correspond to Ly« (top), Ly/3 {middle), and Ha {bottom). 

Notice that there is no significant difference below 1800 km where Lya is isotropic and the 2^3/2 level is unpolarized. 
Above 1800 km the only upper level that is seriously affected is 8^5/2, due to the strong 2^3/2-3^5/2 optical transition. 
The level 8^3/2 cannot be affected directly by the presence of atomic polarization in the 2^3/2 level because the 
transition 2p3/2-8p3/2 is forbidden. However, the fractional alignment of this 8733/2 level is slightly larger when the 
atomic polarization of the lower-level is taken into account, because of the strong coUisional coupling with the 8(^5/2 
and 3(^3/2 levels. The level 8^3/2 is not strongly radiatively coupled to 2^3/2, not even coUisionally to the 8^3/2 level, 
and it is thus not noticeably affected by the presence of atomic polarization in the 2^3/2 level. 

In conclusion, since the atomic polarization of the 2^3/2 level only affects the atomic polarization of the upper 
levels of Ha in the uppermost layers of the model atmosphere, where Ha is already optically thin, the emergent Q/I 
profile of Ha has nothing to do with lower-level polarization. The same applies to the real solar chromosphere (see 
IStepan fc Truiillo Buencill2010[ ). In different scattering environments it might be possible that the Ha scattering polar- 
ization profile is significantly altered by lower-level polarization. This interesting problem will be carefully addressed 
in a forthcoming investigation. 

4.5. The role of collisions 

In semi-empirical m odels of the solar chromosphere (e.g., IVernazza et al.l 119811 : iFontenla et al.l Il998l 120071 : 
lAvrett fc Loeseil 120081 ) the electron and proton densities typically vary between 10^° and IG^^cm"'^ in the upper 
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Fig. 8. — Qjl amplitudes for a close to the limb LOS (/i = 0.1) in the non-magnetic model atmosphere versus the perturbers' density 
Afpcrt, assuming Ne = Np = Npcrt- 

chromosphere, where the Une-center features of the hydrogen hnes studied here originate. In this section, we discuss 
the influence of colhsions with electrons and protons on the emergent / and Q/I profiles. More information on the 
inelastic and depolarizing coUisional rates of our calculations can be found in Appendix [Bl 

We begin by investigating the impact of coUisional transfer between the n — S sublevels on the fractional alignment 
of the hydrogen atomic levels. Obviously, if collisions among the sublevels of n = 3 are neglected then the transfer of 
population and alignment between the ensuing fine structure levels is no longer possible. The top right panel of Fig. [S] 
shows the results of this numerical experiment. Note that when collisions among the sublevels of n = 3 are neglected 
the 3dj levels are not coupled to the 3pj levels. In contrast to the models in which the An = collisions are taken 
into account, the fractional alignment of the 3dj levels increases in the intermediate layers of the model atmosphere 
where Ha starts to become optically thin but the Lyman lines are still opaque. 

The bottom right panel of Fig. |6] helps to clarify the effect of the 2^3/2 polarization on that of the level. If 

the An = collisions are neglected in n = 3 and the polarization of the 2p3/2 level is artificially suppressed, the 
polarization of the 3^5/2 level in the uppermost atmospheric layers is no longer affected by the lower level polarization 
(compare with the other panels of the same figure). The polarization of the 3^5/2 and 3^3/2 levels is due to the 
anisotropy and polarization of the Ha radiation and it is practically uncorrelated with the anisotropy and polarization 
of the Ly a and Ly /3 lines. 

We now turn to showing in Fig. [7] the emergent line intensity and the fractional linear polarization for increasing 
values of the perturbers density. Since the model atmosphere is isothermal, the line source functions mostly decrease 
with height in the atmosphere (see Fig. 14]) and the emergent intensity profiles appear in absorption. The anisotropy 
of the spectral line radiation here is dominated by the limb darkening of the outgoing radiation and the emergent Q /I 
signals are positive (i.e., parallel to the surface of the atmosphere). 

Fig. [7] shows that the sensitivity of the Q/I profile of the Ly a line to coUisional depolarization is relatively small 
(see also Appendix IB. 2p . For densities lower than 10^^ cm""^ there is practically no modification of the emergent Q/I 
profile. A noticeable decrease in the Q/I line-center amplitude can be seen around 10^^ cm^'^, i.e., at a density which 
is considered to be too high for the uppermost layers of the solar chromosphere. 

As shown in Fig. [5] the anisotropy of the spectral line radiation at the atmospheric height where T(/i — 0.1) = 1 
is larger for Ly/3 than for Lya. However, whether or not the amplitude of the scattering polarization Q/I profile 
is larger in Ly/3 than in Lya depends on the perturber's density. The reson is that the Ly/3 line is very sensitive 
to coUisional depolarization because coUisional shuffling of the optical electrons among the fine-structure sublevels is 
much more efficient in n = 3 then in n = 2. Thus, for Alport = lO^^cm"^ the amplitude of the Ly/3 polarization is 
slightly larger than for Lya (see Fig. [7]), while for A'pert = 10^^ cm~'^ is much smaller. 

The most sensitive spectral line to such collisions is Ha (see the right bottom panel of Fig. [7]) . The rate of coUisional 
transitions among the Mj levels becomes comparable to the inverse radiative lifetime of the 3dj levels at around a 
perturbers density of a few times 10^°cm~'^. Note that the source function and anisotropy of the Ha line are less 
sensitive to the thermal structure of the chromosphere than Lya and Ly/3 (because Ha is photo-ionization dominated), 
and that the Ha scattering polarization can be significantly reduced by depolarizing collisions with protons if their 
density exceeds approximately 10^*^ cm~^. 

The previous results are summarized in Fig. [8l which shows the Q/I amplitudes of the three hydrogen lines as a 
function of electron and proton density. 
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Fig. 9. — Left panels: CLV of the emergent Q/I profiles of Ly« (top), Ly/3 (middle), and Ha {bottom). The LOS direction cosine, 
l_i = cos 6, varies from 0.1 (solid lines) to 1 (dash-triple dotted lines). Right panels: CLV of the ensuing line-center polarization amplitudes. 

4.6. Center-to-limb variation ( CLV) of the line profiles 

The left panels of Fig. [9] show the fractional linear polarization profiles calculated for various line-of-sight inclinations 
using our nominal value for the perturber's density (i.e., Alport = 4 x 10^" cm~'^), while the right panels show the center 
to limb variation of the line-center amplitudes. The scattering polarization of the Ly a line is produced in the uppermost 
layers of the model atmosphere by the atomic alignment of the 2p3/2 level that is induced by anisotropic radiation 
pumping in the Lya transition itself (i.e., lsi/2-2p3/2). The amplitude of the emergent Q/I signal can be estimated 
using Eq. [3]from the fractional alignment pg(2p3/2)/po(2p3/2) at the height where the line-center optical depth along 
the LOS is unity (cf. the top left panel of Fig. |6]). Given that the fractional alignment increases with height, the 
maximum polarization amplitude is found close to the limb (see the top right panel of Fig.[5|). 

The formation of the Ly /3 line is affected by its coupling with Ha through the n = 3 levels. This affects mainly the 
wings of the Ly/3 line which show an increase of the Q/I signal (see the middle left panel of Fig. |9]). The wings are 
formed in the region where the line center optical depth of Ly f3 is much higher than unity, but the polarization of the 
level 3p3/2 is non-zero (see the discussion in ^^1^. This effect is more obvious for a LOS with a large n value, i.e., 
closer to the disk center, where the ratio of the Q/I signal in the line center and in the wings becomes smaller due 
to a relatively higher contribution of the deeper layers between 1100 km and 1600 km. Since we only consider thermal 
Doppler broadening in our model, the resulting profiles are quite narrow. If the natural and Stark broadening were 
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Fig. 10. — Fractional alignment of the fine-structure levels in a uniformly magnetized model atmosphere. The perturber's density is 
= 4 X lO^^cm"^. Le/t: Horizontal magnetic field with random azimuth and strength _B = lOG. Right: Horizontal magnetic field 



with random azimuth and strength B = 100 G. 




Fig. 11. — The emergent Q/I profiles of hya (left), Ly/3 {middle), and Ho (right) calculated for a LOS with fi = 0.1. The magnetic 
field is uniform, with a random azimuth XB and a constant inclination 9b = 90°. Three magnetic field strengths are considered: OG {solid 
lines), 10 G {dashed lines), 30 G {dash-dotted lines), and 100 G {dotted lines). 

taken into account, one would obtain broader Q/I wings in Ly/3. We note, however, that even though the partial 
redistribution effects in Ly/3 are weaker than those of Lya (e.g., V ernazza et aLulQSltiHubenv fc Lites.l995i) . coherent 
scattering in the wings could significantly modify the resulting line profiles. 

5. THE MAGNETIZED CASE 

The modification of the emergent linear polarizatio n by the H anle effect provides a tool for dia gnostics of solar and 
stellar magnetic fields (Stenflo 1994; La ndi DegPInnocenti fc Lan dolfl 20041 : iTruiillo Buenoll2001| ). In this section, we 
discuss the effects of magnetic fields, both deterministic and with random azimuth, on the scattering polarization of 
hydrogen lines. 

5.1. Micro-structured magnetic field 

Here we consider the micro-structured field case in which the intensity B and the inclination 9b are fixed at each 
height of the atmosphere but the azimuth xb is random with a uniform distribution in [0, 27r) at scales smaller than the 
mean free path of the line-center photons. The fact that under such circumstances the magnetic field is cylindrically 
symmetric with respect to the vertical axis implies that the rotational symmetry of the problem is preserved, so that 
the Pq^o coherences vanish throughout the atmosphere. The problem is thus numerically similar to the case of a 
non-magnetic atmosphere (see Appendix [C] for details) . The difference is that the micro-structured field decreases the 
alignment of the levels. The Q/I profiles are thus depolarized and the Stokes-J7 parameter remains identically zero. 

It is sufficient to consider two magnetic field strengths, namely 10 and 100 G. It follows from Table [T] that a lOG 
field is sufficient to modify the alignment of the 3p3/2 level and especially those of the 8^3/2 and 8^5/2 levels. The 
2p3/2 level should be only slightly modified because its critical Hanle field, Bh, is about 53 G (i.e., well above 10 G). 
A 100 G field will significantly alter the atomic polarization of the 2^3/2 level while the 3ij levels should already be 
in the Hanle effect saturation regime. Fig. (TOj shows the variation with height in our stellar atmosphere model of the 
fractional alignment of the hydrogen levels, for the case of a random-azimuth horizontal magnetic field with such field 
strengths. 

The left panel of Fig. (TOj shows the results for the 10 G case, using the same scale as in Fig. \6\ As expected, 
Pq(2p3/2)/Po (2^3/2) is decreased but since the critical field of this level is five-times larger than the applied mag- 
netic field, the change is small. The strongest depolarization occurs for the 3dj levels whose alignment is decreased 
by about a factor two with respect to the non-magnetic model. At heights lower than about 1800 km we have 
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From top to bottom: Lyo, Ly/3, and Ha at the indicated /^-value of the LOS. The positive Q-direction is the perpendicular to the stellar 
radius through the observed point. Note that for the forward scattering case the positive Q-direction is the perpendicular to the horizontal 
magnetic field. 

|po(3d5/2)/Po(3(i5/2)| < |/5o(3d3/2)/Po(3d3/2)|, i-e., the fractional polarization of the 8(^3/2 level exceeds that of 3^5/2- 
The reason is that the 8^5/2 level is sensitive to weaker magnetic fields than the 3^3/2 level (see Table [ij. The dashed 
hnes of Fig. [TT] show the emergent Qjl profiles corresponding to the 10 G case. As expected, the Qjl profile of Lya 
is similar to the non-magnetic one (see the solid lines). On the other hand, the Ly/3 and Ha lines become noticeably 
depolarized. 

The right panel of Fig. [TU] considers the 100 G case. Note that now the levels are significantly more depolarized, to 
the extent that the Hanle effect in the n = 3 levels becomes virtually saturated. The alignment of the 2^3/2 level (i.e., 
the only Lya level that can be aligned) is also significantly reduced. The fractional alignment of the Sdj level and of 
the 3p3/2 level is approximately a factor four smaller than in the non- magnetic model (as expected for the saturation 

field regime, see iManso Sainz et al.l I2006I ) and the relation |po(3c?5/2)/Po(3d5/2)| > |Po(3'^3/2)/Po(3c^3/2)| is satisfied 
again throughout the atmosphere because the field intensity is much larger than the critical fields of both of Zdj levels. 
The emergent Qjl profiles of the 100 G model are given by the dotted lines in Fig. [TT] As expected, the amplitudes 
of the Ha and Ly /3 signals are a factor four smaller than in the non-magnetic model. 



5.2. Deterministic magnetic field 
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Fig. 13. — The emergent Qjl profiles of Ly « (left), Ly/3 {middle), and Ho (right) calculated at the disk center (fi = 1). The magnetic 
field in the atmosphere is uniform with the fixed azimuth xb = 90° ^^nd the constant inclination dg = 90°. Three magnetic field strengths 
are considered: 10 G {dashed lines), 30 G {dashed- dotted lines), 100 G {dotted lines), and 200 G {solid lines). 

The positive reference direction we have chosen for Stokes Q is the perpendicular to the stellar radius through the 
observed point. Therefore, it is clear that the Stokes U parameter is zero if the magnetic field has a random-azimuth 
distribution, but non-zero in general if the atmosphere is permeated by a deterministic magnetic field vector (i.e., with 
a well-defined inclination and azimuth). This is illustrated in Fig. I12[ which shows the center to limb variation of the 
Ql I and U /I profiles for the case of a horizontal magnetic field of 100 G with azimuth xb = 0° (i.e., pointing towards 
the observer for a LOS with /i = and perpendicular to the LOS for a /i = 1 disk center observation). 

Note in Fig. [T2]that for the particular case of forward scattering geometry (i.e., LOS with ^ = 1) we have U/I — 0, 
while the Q/I signal is the largest. The fact that for this case U/I = can be easily understood by symmetry reasons. 
Fig. [T51 illustrates how the forward scattering polarization signals change with the field strength. As expected, the Ha 
forward scattering signal is rather small, due to the low radiation anisotropy around the atmospheric height where the 
line center optical depth is unity. On the other hand, the fractional linear polarization of Ly /3 is very significant and 
the polarization profile is quite broad due to the significant alignment of the 3p3/2 level between 1200 and 1700 km. 
In the limit of strong magnetic field, the polarization of Ly a is about a factor two smaller than that of Ly /? and, in 
contrast to Ha and Ly /3 it appears at higher field strengths. 

5.3. The impact of magnetic field gradients 

Whereas the Lya and Ly/3 Stokes Q signals result from just one polarized transition (lsi/2~2p3/2 and lsi/2- 
3p3/2, respectively, cf. Fig. [2|), the scattering polarization of Ha is the result of the superposition of 5 fine-structure 
transitions between the n — 2 and n = 3 j-levels that can, in principle, contribute to the line polarization. As we 
have seen, four of these transitions (i.e., the ones whose upper levels can be aligned) play the dominant role on the 
Ha scattering polarization. In solar- like atmospheres the contribution of the dichroic-only transition 2p3/2-3si/2 can 
be ignored because the atomic polarization of the 2p3/2 level is negligible in the Ha formation region (see ^4.41) . The 
wavelength separation of the Ha line components is much larger than those of the Lyman lines (see Fig. [2]). In the 
presence of magnetic field gradients the non-negligible wavelength separation of the Ha components gives rise to a 
wavelength-dependent Ha nle effect, which is of diagnostic interest because it may produce asymmetries in the emergent 
linear polarization profile (jStepan fc Truiillo Buen"oll2010[ ). In this subsection, we emphasize further the role of these 
overlapping line components on the radiation transfer problem of scattering polarization. 

The individual Ha components are displaced from the line center because of the fine-structue splitting and the Lamb 
shift of the levels. A typical separation of the Ha components is 0.1 A (see Fig. [5]). The thermal Doppler width of 
the line (which is typically about 0.3 A in the solar chromosphere) is significantly larger and, consequently, the Ha 
intensity profile resembles the Gaussian bell. Since the fine-stucture transitions i-j of the Ha line are shifted from 
line center, the optical depth scale resulting from the opacity due to all possible absorptions at a given frequency is 
not symmetric with respect to the central wavelength of any given component, Xij. Rougly speaking, the photons at 
wavelengths — AA and Xij + AA originate from different geometrical depths in the atmosphere. Neglecting selective 
absorption (because Pq (2^3/2) ~ 0), the Stokes Q signal of the emergent radiation along a given LOS can be expressed 
as 

/>oo 

Q(A)- / dtSQ(\,t)e-' , (13) 
Jo 

where t is the optical path along the LOS at the wavelength A, and Sq — sq/tii is the source function of Stokes Q. 
The emergent Stokes Q-profile can be decomposed as 



Q(A)=/ di^^g(A,Oe-* = ^Qy(A), 



(14) 



where Qij and Sq = f^Q/Tji ^re the profile components of the individual transitions and their source functions, 
respectively. We emphasize that 77/ is the total absorption coefficient resulting from all the Ha transitions, whereas 
eg is the emissivity of the i-j transition only. 
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Fig. 14. — The Q/I profiles of H« and of its four main polarizing components calculated for a LOS with ^ = 0.1 and the influence of a 
magnetic field. Top left: The case of a non-magnetic atmosphere. Top right: The case of a micro-structured magnetic field with B = 10 G 
and inclination 9b = 65°. The two bottom panels illustrate the effect of gradients in the magnetic field intensity. Bottom left: Variation 
of the magnetic field strength and fractional alignment of the levels against height in the model atmosphere. The chosen magnetic field 
strength increases with height (Z) according to the sigmoid function B(Z) = Bo/[l + exp{—{Z — Zq) /w)], where Bg = 50 G, Zq = 1300km, 
and w = 100 km. Bottom right: Q/I profile of the emergent Ha line. 

Let us now consider a particular point on the LOS parameterized by the geometrical path s. The emission coefficient 
Cq is symmetric with respect to the component central wavelength Ay , i.e., eg(Ay — AA, s) = eg (Ay -I- AA, s). From 
now on, we will denote this quantity by eg (Ay- ± AA, s). The total absorption coefficient 77/ is only approximately 
symmetric with respect to the line-center wavelength. On the other hand, it is not symmetric with respect to the Ay- 
wavelengths, i.e., ?7/(Ay — AA, s) ^ r\i(\ij + AA, s). Changing the integration variable from t to s, it follows that the 
difference 



,(Ay - AA) 



y(Ay +AA) 

- /q ds T]i{\ij-A\,s') 



ds eg (Ay- ± AA, s) 

Jo ds' T,iiXi,+A\,s')' 



(15) 



is generally non-zero and that the component Qy of the line is asymmetric with respect to its central wavelength Ay . 
We may thus conclude that (1) the Qij/I profiles are asymmetric even in the absence of magnetic field and (2) the 
formation depth of the Qij/I components is a non-trivial function of wavelength. 

The top left panel of Fig. [T4l shows the Qij/I components and the full Q/I profile calculated for a LOS with /i = 0.1 
in our isothermal model atmosphere, assuming B = gauss. A comparison with the transitions central wavelengths 
(see right panel of Fig. [5]) shows that the positions of the Qij/I maxima are shifted from the line center towards the 

wings. This fact can be understood by considering the component's source function Sq{X) = eQ{X)/rii{X) at any point 
in the atmosphere. One realizes that the maximum of Sq (A) is always shifted towards the wing with respect to the 
maximum of eg(A) because the denominator ?7/(A) decreases towards the wing. Note that in the absence of magnetic 
fields the Q/I profile of Ha is clearly dominated by transitions "1" and "2". 

The top right panel of Fig. [14] shows the results of a similar calculation but taking into account the Hanle effect 
of a micro-structured magnetic field of 10 G with an inclination 6b = 65°. The Q/I profile is now skewed towards 
the blue part of the spectrum. This effect can be understood by considering the Qij/I components. In contrast to 
the non-magnetic case (see left panel of Fig. [T4|) . transition "2" now dominates the Q/l profile because the alignment 
of the 3(i5/2 level is strongly depolarized by the Hanle effect. Component "2" contributes mainly in the blue part of 
the line. The contribution from transition "1" in the red part of the profile is smaller and transition "3" decreases 
the Q/I signal (note that transitions "2" and "3" have the same upper level 3^3/2, which is less depolarized than the 
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upper level of transition "1"). Consequently, the whole profile is slightly skewed towards the blue wavelengths. For 
increasingly stronger fields, the importance of the difference between the critical field values of 8^3/2 and 8(^5/2 becomes 
increasingly negligible and the relative Qij/I contributions corresponding to the non- magnetic model are eventually 
restored for field strengths larger than the Ha Hanle saturation field. 

From the previous discussion related to the top panels of Fig. [14] we point out that neither the emergent Q profile 
nor Ql I show any significant asymmetry, even though the individual components are asymmetric. The reason is that 
the shapes of the four Qij/I components of Ha are smooth "skewed gaussians" located aroimd the line center, whose 
superposition gives rise to a rather smooth Q / 1 profile without any noteworthy asymmetry. 

However, something interesting happens if we have magnetic field gradients in the line-core formation region of the 
Ha line (see the bottom panels of Fig.ll4p. The above-mentioned second conclusion provides the clue for understanding 
why the emergent Q / 1 profile shows now a line core asymmetry (LCA). In agreement with our numerical experiments, 
the LCA can be created by a modification of the Qij/I profiles of some components (mainly the "1" and "2" transitions) 
induced by a spatially varying magnetic field. The superposition of the blended linear polarization profiles can then 
give rise to a sizable LCA whose shape depends sensiti vely on the spatia l distribution of the chromospheric magnetic 
field. The interpretation of the Q/I profile observed bv lGandorfeil ()20QCl[ ) in the Ha line led us to suggest that there is 
a significant and abrupt magnetization in the upper chromosphere of the quiet Sun (jStepan fc Truiillo Buen"oll2010l ). 

6. CONCLUDING COMMENTS 

In this first paper of a series on the scattering polarization of hydrogen lines in weakly magnetized stellar atmospheres 
we have considered the case of an exponentially stratified isothermal model atmosphere, which has allowed us to gain 
physical insight on the atomic level polarization produced by optical pumping processes and the operation of the Hanle 
effect in the Lya, Ly/3 and Ha lines. In order to be able to investigate this problem taking into account multilevel 
radiative transfer effects, the overlapping of the fine-structure transitions contributing to each spectral line and the 
impact of collisions and the Hanle effect, we have applied the efficient numerical methods outlined in Appendix \K\ and 
the method explained in Appendix ICl when dealing with micro-structured magnetic fields. 

It may be useful to emphasize the following points: 

• The fractional scattering polarization of Lya, which is solely due to the selective emission of polarization com- 
ponents resulting from the atomic alignment of the 2^3/2 level, can be safely calculated by solving the radiative 
transfer problem assuming a three-level atomic model, with the lsi/2 ground level and the 2pi/2 a-nd 2^3/2 upper 
levels, taking into account the overlapping between the two possible radiative transitions and neglecting collisions 
between the n — 2 fine-structure levels (if the perturber's density iVpert^lO"'^"'^ cm~^). Via the Hanle effect the 
scattering polarization of the Ly a line is sensitive to inclined magnetic fields with strengths 10 G < B < 250 G, 
approximately. 

• Reliable calculations of the fractional scattering polarization of Ly /3, which is solely due to the selective emission 
of polarization components resulting from the atomic alignment of the 8^3/2 level, require taking into account 
the fine-structure levels of the first three hydrogen n-levels, including the collisional transfer of aligment between 
the n = 8 levels. The Hanle effect sensitivity of Ly/3 lies between 8G and 80 G, approximately. 

• Modeling the scattering polarization of Ha, which is in general the result of the contributions from five polarizing 
blended transitions, requires taking into account also the fine-structure levels of the first three hydrogen n-levels, 
including the collisional transfer of aligment between the rt = 8 levels. In solar-like atmospheres the emergent 
linear polarization has nothing to do with lower-level polarization, because in the atmospheric region where 
the Ha polarization is produced the atomic alignment of the 2^3/2 is negligible. Interestingly, the wavelength 
separation between the overlapping Ha transitions makes possible a wavelength-dependent Hanle effect which 
may give rise to asymmetric Q/I profiles when in the presence of magnetic field gradients. The scattering 
polarization of Ha reacts to fields between 1 G and 50 G, approximately. 

In a forthcoming Letter and in the next paper of this series we will describe in great detail the results of our 
calculations in semi-empirical models of the solar atmosphere, with predictions on the scattering polarization signals 
produced by atomic level polarization and the Hanle effect. 

We are grateful to Luca Belluzzi (lAC) and Egidio Landi Degl'Innocenti (University of Firenze) for carefully re- 
viewing the paper and for suggesting various useful improvements. Financial support by the Spanish Ministry of 
Science and Innovation through project AYA2010-18029 (Solar Magnetism and Astrophysical Spectropolarimetry) 
and CONSOLIDER INGENIO CSD2009-00088 (Molecular Astrophysics: The Herschel and Alma Era) is gratefully 
acknowledged. 

APPENDIX 

A. THE NUMERICAL METHOD OF SOLUTION 

In order to obtain the self-consistent solution of the hydrogen scattering polari zation proble m taking into account 
the Hanle effect in multilevel atomic models we have applied a computer program ()Stepanll2008^ based on an accurate 
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Fig. 15. — Relative change versus iteration number it defined as \pQ{it) — PQ{it — l)\/\pQ (it)\, where |pq (it)| = IPq (*' 

euchdian 2-norm of the vector formed by the density matrix multipoles at all spatial grid points i in the unmagnetized model atmosphere. 



The solid lines correspond to the pj] quantities, while the dashed lines to pg. 
structure levels pertaining to the same level n are very similar. 



Note that the relative changes corresponding to the fine 



formal s olver of the radiative transfer equation (iTruiillo Buenol[2003f ) and on a generalization of a very efficient iterative 
scheme ([Truiillo Buenol[l999D to the case of overlapping transitions. 
The Stokes-vector transfer equation (see Eq. ^ can be also written as follows: 



dr 



(Al) 



where dr = —i]jds and the effective source- function vector S'cff = S — K'l being K' — K /rjj — 1 (with 1 the unit 
matrix and S = e/rij). Therefore, the formal solution of the Stokes- vector transfer equation can be expressed as 

I = ^[S K'I]+T , (A2) 

where ^ denotes an operator which reduces to that introduced bv iRvbicki fc Humme"i1 (1992[ ) when atomic level 
polarization is neglected. In contrast to the standard A operator (jRvbicki fc Hummed Il99i ). ^ acts on the emission 
coefficients of a line rather than on its source function. This enables us to treat more easily the present complex 
case of hydrogen spectral lines composed of multiple overlapping transitions. In Eq. (IA2p . we have used the modified 
propagation matrix, K' = K/rjj — 1, and the Stokes source function vector S = e/tji- The coefficient 77/ is the 
standard absorption coefficient of the radiative transfer theory, which resides on the diagonal of K. The vector T 
denotes the contribution resulting from the Stokes vector illuminating the boundary of the integration domain. 
Since ^ is a linear integral operator, we can formally rewrite Eq. (jA2p in the form 



(A3) 



and following the ideas of operator splitting, we wr ite ^ = ^* -I- (^ — >&*) and we use the diagonal of the exact 
operator for the approximate operator yp* (c.f.. IRvbicki fc Hummed Il992[ ). We label the "old" quantities (i.e., 
those calculated using the density matrix elements calculated in the previous iteration) by the superscript and 
we replace the first term in the right-hand side of Eq. (|A3p by 



*[S']^*[S'+] + **[S'-5' 



+1 



(A4) 



The quantities in the remaining terms are taken from the previous iteration, hence we can rewrite the iterative scheme 
as 

I = 1+ + '^*[S - S+]. (A5) 
The source function vectors in the current and in the previous iteration are defined as 

S = -^, S+ = ^. (A6) 
Vi Vi 

Note that the ^-operator only acts linear ly on the emission co e fficient and that the highly non-linear d ependence on the 
absorption coefficient is suppressed (see ITruiillo Buenolll999t iManso Sainz fc Truiillo Buenci|[2"003al , for an analogous 
strategy). In this way the formal solution of the Stokes-vector transfer equation at a given frequency can be expressed 
as a linear function of the density matrix elements of the upper levels of the active transitions. 
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The next step is to use Eq. (jA5[) in the calculation of the radiation field tensors of Eq. pT|) and to use the resulting 
expressions in the statistical equilibrium equations (SEE) . The source function vectors S and S'^ have to be explicitly 
expressed in terms of the density matrix elements of the upper levels of the radiative transitions. For the sake 
of simplicity, we denote in the following the multipole components of the density matrix by pk- The population 
multipoles pjj are the convergence-driving quantities. We will identify them by a bar over the index, for example pm- 
Using only the terms related to the Jg tensor,^ the preconditioning scheme of SEE is of the form 

^0 («J>fc JQ{ij)^Pk + UjjniPmPk - pirPk) ■ (A7) 

Tn 

The summation runs over all the population multipoles pm of the levels emitting at the frequencies of the i-j transition 
under consideration. The pk elements correspond to the lower-level multipole absorbing at the frequencies of that 
transition. The coefficients L^jj^ can be straightforwardly calculated from the a pproximate ^* operator while doing 
the formal solution via the DELOPAR method proposed bv lTruiillo Buenol ()2003() . The advantage of using the diagonal 
approximate operator is that the SEE remain local and it is easy to evaluate Lijm from the interpolation coefficients 
of the DELOPAR formal solution method of the transfer equation (see lMaiiso" "Sah^ ([2OOI ). and lSteoanl (pOOl for the 
explicit form of the Lijjn coefficients of the multi- level atom). 

Whereas the second term in the bracket in the right-hand side of Eq. (jA7[) is linear in the "new" density matrix 
element, the first product is no t. In order to preserv e linearity of SEE, the term Pm.Pk has t o be replaced by Pmpt (see 
the analogous linearization by iTruiillo Buenol Il999t iManso Sainz fc Truiillo Buenol l2003al) . Finally, we end up with 
the iterative scheme 

Jo{ij)Pk Joi^jVPk + X! LijTrziPfnPt - ptrPk) , (A8) 

Tn 

^Q(u>fc^^Q(u)+Pfc, for K>0. (A9) 

As an example, consider the Lya line which consists of two overlapping transitions: 2px/2~l'Si/2 ^md 2p3/2^1si/2- 
The product Jq (2^3/2 — lsi/2)/3o(lsi/2) will thus become 

>/o(2P3/2 - l'Sl/2)^Po(lSl/2) + ■^2p3/2-lsi/2,2!3i/2 
X [pO(2pi/2)p^(lSl/2)+ - pS!(2pi/2) + pJ(lSi/2)] 

+ ^2p3/2-lSl/2, 2p3/2 

X [pS](2p3/2)4(lSi/2)+ - pS](2p3/2) + p'(lSl/2)] . (AlO) 

The resulting SEE, linear in all pj, can be solved using the standard methods of numerical linear algebra. Figure [T5] 
illustrates the convergence behavior of the iterative method for the density matrix elements corresponding to the 
hydrogen model atom of Fig. 1, considering radiative transfer effects in an isothermal model atmosphere. 

In the c ase of negligible pola r izatio n, the previously described technique is equivalent to the full-preconditioning 
strategy of iRvbicki fc Hummeil ( 1992f). For multilevel atomic models w ithout overlapping transitions, the method is 



equivalent to the one applied bv lManso Sainz fc Truiillo Buenol ()2003af) . We note that an analogous iterative scheme 
can be used for modeling atomic lines with even more complicate d structure, such a s those involving levels in the 
incomplete Paschen-Back regime or levels with hyperfinc structure (|Stepanll2009l l2008[ l . 

B. COLLISIONS 

We treat collisions within the f ramework of the irreducible components of the density matrix (e.g., 
iLandi Degl'Innocenti fc Landolfill2004[ ). The collisional rates of transfer and relaxation of the components of the 
atomic density matrix due to interactions with the ambient thermal particles have the form C'^^\n' i' j' ,nlj). In this 
paper, we do not explicitly distinguish between inelastic (i.e., exciting) and superelastic (i.e., de-exciting collisions). 
The C^''^ rates are responsible for transfer and relaxation of the level populations and they are equal to the familiar 
collisional rates of the unpolarized theory, which will be denoted by Cnij^n'tj' in the following. The C*^^^ rates are 
responsible for transfer of alignment between the atomic levels. The effect of elastic depolarizing collisions, quantified 
by the D'^^^ rate (see Eqs. 7.101 in ILandi Degl'Innocenti fc Landolfill2004l ) . is negligible for the case of hydrogen levels 
at the densities and temperatures typical of the upper solar chromosphere. The depolarization of the atomic levels 
is almost exclusively due to the "weakly inelastic" transitions nlj — n(^ ± between the close- lying fine structure 
levels of a given level n, which are induced mainly by collisions with the ambient protons. 

B.l. Inelastic collisions 

A careful consideration of the collisional processes is needed in the unpolarized non-LTE problem and even more in 
the present case of scattering polarization. 

The role of inelastic collisions with thermal electrons is indeed similar in both kinds of models. Unfortunately, 
the inelastic excitation cross-sections a{n n',E) by electrons are typically poorly known. Hydrogen cross-sections 

^ These are responsible for the convergence rate. 
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TABLE 2 
Excitation and de-excitation 

COLLISIONAL RATES [s""' cm"''] FOR 

T = 10* K, ric = 4 X 10"cm-3. Data 

OBTAINED FROM 'P rZYBILLA Hi BuTLe"rI 



1 2 1.01(-2) 3.49(+2) 

1 3 3.62(-4) 5.01(+1) 

2 3 1.63(+3) 6.48(+3) 



Note. — a{6): a x lO''. 

are notoriously difficult to calculate because of the existence of high number of levels and resonances. From author 
to author, the uncertainty of cross-sections and, consequently, the coUisional rates, can vary by tens of percent (see 
iPrzvbilla fc But ler 2004, and references therein). 

If one considers the fine structure of the levels, the situation with the cross-section data is even more uncertain. 
To our knowledge, only the excitation cross- sections for terms fro m the ground level, a{ls — ^ nl,E), are currently 
available from laboratory measurements fe.g.. lJanev Sz Smithlll993[ ). At low energies, which play the most important 
role in the solar atmosphere, these transitions do not obey the strong selection rules of the dipolar optical transitions, 
|A^| = 1, and the non-dipolar excitations such as Is — > 2s, Is — >■ 3s, Is — ^ 3(i, have rates comparable to the dipolar 
ones. No experimental data for the cross-sections between the fine structure levels of the subordinate lines are available 
and we have to rely on theoretical calculations. 

The approach we hav e chosen in this work is the following. We have adopted the most reliable data of 
IPrzvbilla fc Butleil (|2004f ) available nowadays for the excitation rates, C„,^„^, in a wide range of temperatures and 
we have used the approximation that the coUisional rates between the Zeeman sublevels do not depend on the or- 
bital angular momentum, nor on the magnetic quantum number of the ni — >■ transition under consideration. The 
normalization condition implies the following equation for the population, Nn^^, of the upper Bohr level,'' 

^ - E E - E f E ' (Bl) 

ji iu V ji J 



It is easy to show that the excitation coUisional rate between the two fine-structure levels is 

El: 

3n 



^jl^ju — ^^Cni^n^ , (B2) 



and the relaxation rate 

where gj, (gj^) and gn^ is the statistical weight of the lower (upper) fine structure level and of the full n-level, 
respectively.* Note that the rates (jB2[) and (jB3| connect the fine-structure levels whose populations Nj are related to 
the zero-rank of the density matrix element PQ{nlj) = Nj/^2j + 1 (see i 33.1|) . Reform ulation of the previous expressions 
in the framework of irreducible tensorial operators is straightforward (see Sect. 7.13 of lLandi Degl'Innocenti fc Landolfil 
|2004 for details). The coUisional rates between different n- levels that have been used in this paper are tabulated in 
Table El 

We use the approximation described above in order to guarantee that the coUisional rates of different lines are 
consistent. We have numerically verified that replacing the Cisjj^^n^i^j^ rates by those calculated from the available 
excitation cross-sections a{ni£i — >■ n„£„) of iJanev fc SmithI (jl993[ ) does not modify the Lyman line intensities nor the 
polarization by more than a few percent of their amplitude. Such error is quite acceptable taking into account the 
uncertainty of the available coUisional da ta. It is also worth to me ntion that the Balmer lines are even less afi^ected by 
a modification of the coUisional rates (cf. IPrzvbilla fc Butl"ei1l2004l ) , especially in semi-empirical chromospheric models 
where photoionization dominates over the coUisional transitions. 

B.2. Collisions connecting nearby j-levels 

While the intensity spectrum of the hydrogen lines can be modeled neglecting the fine structure of the levels, the 
fine structure and the effects of depolarizing collisions with different atmospheric species (electrons, ions, etc.) have 
to be taken into account when it comes to modeling the linear polarization produced by atomic level polarization. 

The levels of neutral hydrogen show a high degree of accidental (quasi) degeneracy. The fine-structure levels nij and 
n£'j' pertaining to the same Bohr level n are very close to each other, with a typical separation of 10^^ to 10^^ cV for 

For the sake of simplicity we only explicitly write the rates between two n-levels. 
* It is Qj = 2j -I- 1 and Qn = In?" . We have neglected the energy differences among the different fine structure levels pertaining to the 
same n-level in comparison to the separation between different n-levels. 



22 



Jifi Stepan and Javier Trujillo Bueno 



o 
o 



0) 

C 

o 
+^ 

CO 

c 

cc 




CO 


O IT 

3.5 ■ 


"f 


3.0 - 


o 




CO 


2.5 f 


o 






2.0 - 


a> 


i 


1.5 r 


c 




o 


1.0 r 


CO 


0.5 r 


c 


H 


0.0 : 



proton collisions 
[10"' cm' s'] 




2x10 4x10 6x10 
Temperature [K] 



8x10 



2x10 4x10 6x10 
Temperature [K] 



8x10 



Fig. 16. — CoUisional transition rates of the population transfer per perturber as a function of the plasma temperature for all dipolar 
transitions 2£j 2(1 it l)j' {dashed lines) and 3£j — > 3{£ ± {solid lines). The rates have been calculated for a proton and electron 
density of 10^'' cm""^. The curves show the rates calculated using the semi-cla ssical perturbation method of ISahal-Brechot et al.l l l 19961 ) 
between 3 000 K and 90 000 K, which is based on the impact approximation (see IBommier fc Sahal-BrechotI 1 119911) for a discussion on the 
validity conditions. Left: collisions with electrons. Right: collisions with protons. Note the different scales in the vertical axes. 

n = 3. It can be shown that, at the typical densities of the upper solar chromosphere, coUisional transitions induced 
by electrons and mainly by protons, both having a long-range interaction potential, play a significant role on the 
depolarizatio n process, the coUisional rates due to protons bei ng about one order of magnitude hi gher than those due 
to electron s (Bomr nier et al.iri986bl : ISahal-Br echot et aI1ll996f) . Collisions with neutral hvdrogen ()Ali fc Griem|[l965l : 
iBarklenTet al. 200ij), which play a significant role in depolarizing other lines of the second solar spectrum, can be safely 
neglected in the case of hydrogen lines at the considered densities. Collisions with protons and to a lesser exte nt with 
electrons also play a significant role in broadening of the hydrogen lines via the Stark effect (e.g.. IStehlelll996[) .^ Here 
we only consider the depolarizing role of collisions, with protons and electrons, connecting nearby j-levels pertaining to 
each n-level. In general, the depolarizing effect of collisions with Hen ions should also be taken into account, although 
its significance is expected to be much smaller than that due to collisions with protons (cf., Bommier et al. 1986b). 
All these additional ingredients will be considered in subsequent papers. 

Collisions of hydrogen atoms with perturbers having a Maxwellian velocity distribution tend to produce partial or 
full equilibration of the populations of the Zeeman sublevels, i.e., to a reduction of the atomic level polarization. Since 
the thermal energy of the perturbers (~ 1 eV) is by several orders of magnitude larger than the transition thresholds, 
the most efficient transitions are the dipolar ones, n£j — > n(£ ± Strictly speaking, these transitions are inelastic 
but they actually have a depolarizing effect and they are also responsible for transfer of alignment between the fine 
structure levels. These collisions redu ce the lifetime of the excite d levels. Consequently, the width of the fine structure 
levels increases (see Tables 4 and 5 of ISahal-Brechot et al.lll996l ). If the density of protons becomes too high (roughly 
about 10^'^ cm~^ for the case of n = 3), the levels start to partially overlap and the quantum coherences between them 
should be taken into account in calculati ons of the polarizat ion produced by radiation scattering. Moreover, the impact 
approximation breaks at high densities (jStehle et al.lll983l) . But it is worth mentioning that at such a high density of 
perturbers the atomic polarization of n = 3 level would be practically destroyed (cf. Fig. El) regardless of the theory 
used. In fact, such high densities are only expected in the deep layers of the atmosphere (where also Hi-Hi collisions 
have to be taken into account) and affect mainly the line wings. In the upper chromosphere the fine structure levels 
can be considered as being well separated. 

The dipolar c ross-se ct ions can be calculated using the semi-classica l perturbation theor y (see ISeatonI (|1962t ). 
IBommier et al.l ()1986bt ). IBommier fc Sahal-BrechotI (|1991| l. and mainly ISahal-Brechot et all ()1996l ) for a detailed 
derivation of the coUisional rates in the impact approximation^"). Using this approach, we have calculated the mul- 
tipolar components of the coUisional transition rates (7*^°^ and C*^^^ at each atmospheric height and for every allowed 
transition. The coUisional rates C^*^^ by electrons and protons (normalized to the unit de nsity of the perturbers) 
between the j-levels of n = 2 and 3 calculated using the theory of ISahal-Brechot et al.l (|1996l ) can be found in Fig. \W[ 
Multiplication of these quantities by the actual perturbers density gives the standard coUisional rate in units of s~^. 
When compared to the inverse lifetime of the level, it quantifies the degree of coUisional depolarization of the level. 
As seen in Fig. I16| collisions with protons are about one order of magnitude more efficient than those with electrons. 
It is also clear that the .7'- levels of n = 3 are depolarized significantly more than those of level n = 2 (cf. Table 2 of 
ISahal-Brechot et al.lll996h . This is why the Lya scattering polarization is less affected by the depolarizing collisions 
than Ha and Ly/3 (see right panels of Fig. |7]). Using the data of the right panel of Fig. [TBI reveals the reason why 
coUisional depolarization of Ha becomes critical around a proton density Up « lO^^cm"'^. Using 10~^cm^s~^ as a 
typical value of the coUisional rate per perturber we obtain 10^^ cm^"^ x 10^"^ cm^ s~^ — 10® s^^ which is comparable 
to the inverse radiative lifetime of the upper levels of Ha. 



^ Note that coUisional depolarization and the Stark broadening are intimately related processes f see ISahal-Brechotil2009l . and references 
therein). 

Note that the authors use a notation slightly different from that of lLandi Degl'Innocenti fc Landolfil l|2004l ) . 
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Another important phenomenon related to the depolarizing effect of such collisions connecting nearby j-levels is a 
modification of the critical Hanle field of any given level. Since the value of the critical Hanle field, Bh, increases with 
the inverse lifetime of the level (cf. Eq. [T|) the higher collisional rates make the level less sensitive to the magnetic field. 
As mentioned above, the collisional rates between the j-levels of level n = 3 at a typical chromospheric density can be 
of the order of the inverse radiative lifetime. Consequently, the critical Hanle field can easily be doubled by collisional 
quenching. As a result, the amplitude of the emergent polarization is smaller due to collisional depolarization and the 
onset of the Hanle depolarization is shifted towards higher intensities of the magnetic field. 

C. MULTILEVEL ATOM IN THE PRESENCE OF A MICRO-STRUCTURED MAGNETIC FIELD WITH RANDOM AZIMUTH 

A well-known strategy for modeling the Hanle depolarization of various solar spectral lines is the microturbulent 
magnetic field approximatio n in which the spatial scale of the field var iation is smaller than the mean free path of 
the line-center photons (e.g.. IStenflolll994 iLandi Degl'Innocenti fc: Land olfi 200 4)) . Such a model has been developed 
further by several authors for the case of a two-level atom model ("Truiil lo Bueno fc Mans o Saindll999[ ). for multilevel 
atoms in the Hanle effect saturation regime (Manso Sainz ct al. 2006) and in the optically thin regime (Bclluzzi et al.l 
l2007f) . In this appendix, we consider instead the case of a single- valued micro-structured magnetic field having a fixed 
(but in general height-dependent) inclination 6b and a uniformly distributed azimuth xb at microscopic scales. We 
show how to handle easily such a micro-structured magnetic field when solving the non-LTE problem of the second 
kind for a multilevel atom in a cylindrically symmetric atmosphere with any degree of radiation field anisotropy. 

Let us consider a statistical ensemble of multi-level atoms in a magnetized plasma. We suppose that the magnetic 
field azimuth xb random with a uniform distribution while both the inclination 6b and the intensity B are fixed at 
any given point. Physically, this approximation corresponds to a situation in which the magnetic field varies rapidly 
on a geometrical scale much smaller than the mean free path of the photons. The density mat rix of an en semble of N 
non-interacting atoms is simply an average of the density matrices pi of the individual atoms ()Fandll957t ). 

i 

For the sake of simplicity, we do not explicitly consider the time dependence of the quantities. Every single atom 
is under the infiuence of a deterministic magnetic field. The density matrix of a subset of atoms being embedded 
in the magnetic field with azimuth xb will be denoted by p{xb)- The density matrix of the whole ensemble is a 
superposition of the density matrices of these sub-ensembles weighted over the field distribution. We make a transition 
from discrete to continuous quantities and express the density matrices in the basis of irreducible tensorial operators 
in the atmospheric reference frame with the quantization axis Z being parallel to the vertical axis of the atmosphere. 
Assuming that the distribution of the field's azimuth xb is uniform in the [0, 27r) interval, we can calculate the local 
density matrix 

[pq]a = ^1 dXB [p^iXB)]^, (C2) 

where the subscript A refers to quantities expressed in the atmospheric reference frame. 

Rotation of an "old" reference frame to a "new" one by an angle x around the Z axis leads to the following 
transformation of the density matrix elements (jLandi Degl'Innocenti fc Landolfill2004l ) 

[p^U^^J^iPpUd-D^giR), (C3) 
p 

where VpQ^R) is the rotation matrix corresponding to the rotation of the reference frame parametrized by the Euler 
angles, R — (0,0, x)- In our case, the rotation matrix elements are simply 'Dpq{R) = e~^^'^5pQ] thus 

[p§']ncw = e-*^'? [peloid. (C4) 

Using this transformation in Eq. (|C2[) leads to 

[pq]a -hT'^^^ [PQiXB)]^ = i/^^^^^'^'"'' [PQiXB)]^^ , (C5) 

where [Pq {xb)]xb is expressed in the reference frame in which the magnetic field has zero azimuth. 

Taking into account that any one-dimensional model atmosphere permeated by such a micro-structured field implies 
that the radiation field is cylindrically symmetric with respect to the vertical, it follows that the system is cylindrically 
symmetric as a whole and that [Pq {xb)]xb = Pq{^)j being a solution of the statistical equilibrium equations in the 
case of zero-azimuth magnetic field, is independent of xb- The integral of the exponential in Eq. (jCSP produces the 
Kronecker 6qq. Finally we have 

[pq]a-Pq(^)^Qo- (C6) 

Eq. (jC6[) gives a simple strategy for the numerical solution of the problem. At every point in the atmosphere, 
the statistical equilibrium equations have to be constructed as if the magnetic field was deterministic with azimuth 
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Xb = 0.^^ Solution of the equations gives the Pg (0) density matrix elements. The desired density matrix is then given 

by Eq. (IC6[) . i.e., by the Po (0) elements. These elements can then be used for calculating the local radiative transfer 
coefficients and the RTE can be solved to obtain the emergent fractional polarization profile. 
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In fact, the azimuth can be arbitrary because the pQ components do not depend on it. 



